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Neutron star properties, such as its mass, radius, and moment of inertia, are calculated by solving 
the Tolman-Oppenheimer-Volkov (TOV) equations using the ring-diagram equation of state (EOS) 
obtained from realistic low-momentum NN interactions Viow-k- Several NN potentials (CDBonn, 
Nijmegen, Argonne V18 and BonnA) have been employed to calculate the ring-diagram EOS where 
the particle-particle hole-hole ring diagrams are summed to all orders. The proton fractions for 
different radial regions of a /3-stable neutron star are determined from the chemical potential con- 
ditions — = Me = Mm- The neutron star masses, radii and moments of inertia given by the 
above potentials all tend to be too small compared with the accepted values. Our results are largely 
improved with the inclusion of a Skyrme-type three-body force based on Brown-Rho scalings where 
the in-medium meson masses, particularly those of uj, p and a, are slightly decreased compared 
with their in-vacuum values. Representative results using such medium corrected interactions are 
maximum neutron-star mass 1.8Mq with radius 9 km and moment of inertia ~ 60Afofcm^, 
values given by the four NN potentials being nearly the same. The effects of nuclei-crust EOSs on 
properties of neutron stars are discussed. 
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I. INTRODUCTION 



Neutron stars are very interesting physical systems and 
their properties, such as masses and radii, can be derived 
from the equation of state (EOS) of the nuclear medium 
contained in them. In carrying out such derivation, there 
is, however, a well-known difhculty, namely the EOS is 
not fully known. Determination of the EOS for neutron 
stars is an important yet challenging undertaking. As 
reviewed in [if, 0, H, 0, Hi, this topic has been exten- 
sively studied and much progress has been made. Gener- 
ally speaking, there are two complementary approaches 
to determine the EOS. One is to deduce it from heavy- 
ion collision experiments, and crucial information about 
the EOS has already been obtained [J, H, 0, Hi • Another 
approach is to calculate the EOS microscopically from a 
many- body theory. (See, e.g. 0, i, ^ and references 
quoted therein.) As is well-known, there are a number 
of difhculties in this approach. Before discussing them, 
let us first briefly outline the derivation of neutron-star 
properties from its EOS. One starts from the Tolman- 
Oppenheimer-Volkov (TOV) equations 
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where p{r) is the pressure at radius r and M(r) is the 
gravitational mass inside r. G is the gravitational con- 
stant and e(r) is the energy density inclusive of the rest 
mass density. The solutions of these equations are ob- 
tained by integrating them out from the neutron-star 
center till its edge where p is zero. (Excellent pedagogi- 
cal reviews on neutron stars and TOV equations can be 
found in e.g. 0, [HI.) In solving the above equations, 
an indispensable ingredient is clearly the nuclear matter 
EOS for the energy density e(n), n being the medium 
density. As the density at the neutron star center is typi- 
cally very high (several times higher than normal nuclear 
saturation density of no — 0.16/m~'^), we need to have 
the above EOS over a wide range of densities, from very 
low to very high. 

In the present work we shall calculate the nuclear EOS 
directly from a fundamental nucleon-nucleon (NN) in- 
teraction Vnn and then use it to calculate neutron star 
properties by way of the TOV equations. There have 
been neutron-star calculations using a number of EOSs, 
most of which empirically determined, and the mass- 
radius trajectories given by them are widely different 
from each other (see, e.g. Fig. 2 of [l|). To deter- 
mine the EOS with less uncertainty would certainly be 
desirable. There are a number of different NN potential 
models such as the CDBonn [l2|, Nijmegen [l3|, Argonne 
V18 [3 and BonnA fl^ potentials. These potentials all 
possess strong short-range repulsions and to use them 
in many-body calculations one needs first take care of 
their short-range correlations by way of some renormal- 
ization methods. We shall use in the present work the 
recently developed renormalization group method which 
converts Vn n into an effective low-momentum NN inter- 
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action Viow-k O H [11,120,1211. An advantage of 
this interaction is its near uniqueness, in the sense that 
the Viow-kS derived from different reahstic NN interac- 
tions are nearly the same. Also Viow-k is a smooth po- 
tential suitable for being directly used in many-body cal- 
culations. This Viow-k will then be used to calculate the 
nuclear EOS using a recently developed low-momentum 
ring-diagram approach [22i] , where the particle-particle 
hole-hole (pphh) ring diagrams of the EOS are summed 
to all orders. The above procedures will be discussed in 
more detail later on. 

We shall also study the effects of Brown-Rho (BR) 
scaling [2^, l2 7| . Ip sl. l29j on neutron star properties. As 
discussed in |22| . low- momentum ring diagram calcu- 
lations using two-body NN interactions alone are not 
able to reproduce the empirical properties for symmet- 
ric nuclear matter; the calculated energy per particle 
(Eq/A) and saturation density (ng) are both too high 
compared with the empirical values of Eo/A ~-16 MeV 
and no ~ 0.16 fm^^. A main idea of the BR scaling 
is that the masses of in-medium mesons are generally 
suppressed, because of their interactions with the back- 
ground medium, compared with mesons in free space. As 
a result, the NN interaction in the nuclear medium can 
be significantly different from that in free space, partic- 
ularly at high density. Effects from such medium modi- 
fications have been found to be very helpful in reproduc- 
ing the empirical properties of symmetric nuclear mat- 
ter [22]. Dirac-Brueckner-Hartree-Fock (DBHF) nuclear 
matter calculations have been conducted with and with- 
out BR scaling [H, 12^, . In addition, BR scaling has 
played an essential role in explaining the extremely long 
life time of ^''C /3-decay [s^l- As mentioned earlier, the 
central density of neutron stars is typically rather high, 
^ SriQ or higher, no being the saturation density of nor- 
mal nuclear matter. At such high density, the effect of 
BR scaling should be especially significant. Neutron stars 
may provide an important test for BR scaling. 

In the following, we shall first describe in section II the 
derivation of the low-momentum NN interaction Viow-k, 
on which our ring-diagram EOS will be based. Our 
method for the all order summation of the ring diagrams 
will also be addressed. Previously such summation has 
been carried out for neutron matter [sij and for symmet- 
ric nuclear matter (2^ . In the present work we consider 
/3-stable nuclear matter composed of neutrons, protons, 
electrons, and muons. Thus we need to calculate ring dia- 
grams for asymmetric nuclear matter whose neutron and 
proton fractions are different. An improved treatment for 
the angle averaged proton-neutron Pauli exclusion oper- 
ator will be discussed. In section III we shall outline 
the BR scaling for in- medium NN interactions. We shall 
discuss that the effect of BR scaling can be satisfacto- 
rily simulated by an empirical three-body force of the 
Skyrme type. In section IV, we shall present and discuss 
our results of neutron star calculations based on ring- 
diagram pure-neutron EOSs, ring-diagram /3-stable EOSs 
consisted of neutrons, protons, electrons and muons, and 



the well-known EOS of Baym, Pethick and Sutherland 
(BPS) [32] for the nuclei crust of neutron stars. 



II. RING-DIAGRAM EOS FOR ASYMMETRIC 
NUCLEAR MATTER 

In our calculation of neutron star properties, we shall 
employ a nuclear EOS derived microscopically from real- 
istic NN potentials Vnn- Such microscopic calculations 
would provide a test if it is possible to derive neutron 
star properties starting from an underlying NN inter- 
action. In this section, we shall describe the methods 
for this derivation. A first step in this regard is to de- 
rive an effective low- momentum interaction Viow-k by 
way of a renormalization procedure, the details of which 
have been described in [li, HI Hi HI, HI l2l|. Here 
we shall just briefly outline its main steps. It is gener- 
ally believed that the low-energy properties of physical 
systems can be satisfactorily described by an effective 
theory confined within a low-energy (or low-momentum) 
model space [20j. In addition, the high-momentum (short 
range) parts of various Vnn models are model dependent 
and rather uncertain [20| . 

Motivated by the above two considerations, the fol- 
lowing low-momentum renormalization (or model-space) 
approach has been introduced. Namely one employes a 
low-momentum model space where all particles have mo- 
mentum less than a cut-off scale A. The corresponding 
renormalized effective NN interaction is Viow-k, which 
is obtained by integrating out the fc > A momentum 
components of V/v AT. This "integrating-out" procedure is 
carried out by way of a T-matrix equivalence approach. 
We start from the full-space T-matrix equation 



T{k',k,k^) = VNNik',k) 
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where V denotes the principal-value integration. Notice 
that in the above equation the intermediate state mo- 
mentum q is integrated from to oo. We then define an 
effective low-momentum T-matrix by 

Tion,^k{p\p,P^) = Viow^k{p\p) 

, ^ 2, Viow-k{p' ,q)Tiow-k{q,P,p'^) 
+ V q dq , (3) 

Jo P - q 

where the intermediate state momentum is integrated 
from to A, the momentum space cutoff. The low mo- 
mentum interaction Viow-k is then obtained from the 
above equations by requiring the T-matrix equivalence 
condition to hold, namely 

T{p,p,p^) = Tioni~k{p',P,P^); {p',p) < A. (4) 

Note that T and Ti^^^k are both half on energy shell, and 
they are equivalent within A. The low-energy (< A^) 
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phase shifts of Vnn are preserved by Viow-k and so is 
the deuteron binding energy. As we shall discuss later, 
for neutron star calculations we need to employ a cut-off 
scale of A ~ 3 fm~^. 

For many years, the Brueckner-Hartree-Fock (BHF) 
[S H [la [Hi, III and the DBHF 0] methods were 
the primary framework for nuclear matter calculations. 
(DBHF is a relativistic generalization of BHF). In both 
BHF and DBHF the G-matrix interaction is employed. 
This G-matrix interaction is energy dependent. This en- 
ergy dependence adds complications to calculations. In 
contrast, the above Viow-k is energy independent which 
facilitates the calculation of nuclear EOS. 

We shall use Viow-k to calculate the EOS of asym- 
metric nuclear matter of total density n and asymmetric 
parameter a 

n = nn+np] a = , (5) 

Tlfi 'T' Tip 

where n„ and Up denote respectively the neutron and pro- 
ton density and they are related to the respective Fermi 
momentum by kp^/{3n'^) and kpp/{3TT^). The proton 
fraction is x = (1 — a)/2. 
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(a) (b) (c) 

FIG. 1: Diagrams included in the all-order pphh ring-diagram 
summation for the ground state energy shift of nuclear matter. 
Each dashed line represents a Viom-k vertex. 

In our ring-diagram EOS calculation, the ground-state 
energy shift AE is given by the all-order sum of the pphh 
ring diagrams as illustrated in Fig. 1, where (a), (b) 
and (c) are respectively 1st-, 4th- and 8th-order such 
diagrams. {AEo is defined as (£^0 - E^'''"') where Eq 
is the true ground state energy and E^^'^'^ that for the 
non-interacting system.) Following (2^ . we have 

•^0 m ijkKK 

xY:;^(kl,\){ij\Vio^,^k\kl), (6) 

where the transition amplitudes are given by the RPA 
equation 

^[(ei + <^j)5ijM + KlJ] ~ fifj){ii\yiow-k\kl)] 

kl 

X r„(fc;,A) = w„r„(ij,A); {i,j,k,l)<K. (7) 

In the above, the single particle (s.p.) indices (i, j, ...fc, I) 
denote both protons and neutrons. The s.p. energies e 



are the Hartree-Fock energies given by 

ek^h^k''/2m+ J2 {kh\Vio^^k\kh), (8) 

h<kF{h) 

where kp{h) = kpn if h is neutron and = kpp if it is pro- 
ton. Clearly depends on the nuclear matter density n 
and the asymmetric parameter a. The occupation factors 
fi and fj of Eq.(7) are given by = 1 for k < kp{a) 
and /a = for fc > kpia); also /q = (1 - fa)- Again 
kp(a) = kpn if a is a neutron and = kpp if it is a proton. 
The factor {fifj — fifj) is clearly also dependent on n 
and a. Note that the normalization condition for Ym in 
Eq.(6) is {Ym\^\Yra) = -1 and Q{h, %) = - /,/,) 
(39| . In addition, S in Eq.(6) means we sum over only 

m 

those solutions of the RPA equation ([7]) which are dom- 
inated by hole-hole components as indicated by the nor- 
malization condition. Note that there is a strength pa- 
rameter A in the above, and it is integrated from to 
1. 

The amplitudes Y in Eq.(6) actually represent the 
overlap matrix elements 

F„;(fc/,A) = (*„(A,A-2)|a,afc|*o(A,A)), (9) 

where ^o(A, A) denotes the true ground state of nuclear 
matter which has A particles while 'I'„i(A, A — 2) the 
mth true eigenstate of the (A-2) system. If there is no 
ground-state correlation (i.e. ^'o is a closed Fermi sea), 
we have Y^{kl, A) = fkfi and Eq.(6) reduces to the HF 
result. Clearly our EOS includes effect of ground-state 
correlation generated by the all-order sum of pphh ring 
diagrams. 

A computational aspect may be mentioned. We shall 
solve ring-diagram equations on the relative and center- 
of-mass (RCM) coordinates k and K [i^. (fc = {ki — 
kj)/2 and K = ki + kj.) In so doing, the treatment of the 
Pauh operator Q{ki,kj) = {fifj - fifj) of Eq.(7) plays 
an important role. This operator is defined in the labo- 
ratory frame, with its value being either 1 or -1 (for pphh 
ring diagrams). In our calculation, however, we need Q in 
the RCM coordinates. Angle-average approximations are 
commonly used in nuclear matter calculations, and with 
them we can obtain angle-averaged Q{k,K). For sym- 
metric nuclear matter, detailed expressions for Q have 
been given, see Eqs.(4.9) and (4.9a) of The results 
for this case, where the Fermi surfaces for neutron and 
proton are equal, are already fairly complicated. Deriva- 
tion of the angle-averaged Pauli operator for asymmetric 
nuclear matter has been worked out in [4l[ . Their deriva- 
tion and results are both considerably more complicated 
than the symmetric case. It would be desirable if they can 
be simplified. We have found that this can be attained 
by way of a scale transformation. Namely we introduce 
new s.p. neutron and proton momentum coordinates k'^ 
and k' defined by 



k'n ^ knJkpp/kpn, k'p ^ kpJ kpn / kpp. (10) 
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On this new frame, the Fermi surfaces for neutron 
and proton become equivalent, both being {kpnkppY/'^] 
asymmetric nuclear matter can then be treated effec- 
tively as a symmetric one as far as the Pauli operator 
is concerned. We have found that this transformation 
largely simplifies the derivation and calculation of the 
angle-averaged asymmetric Pauli operator Q and will be 
employed in the present work. 

III. IN-MEDIUM NN INTERACTIONS BASED 
ON BROWN-RHO SCALING 

A main purpose of the present work is to study whether 
neutron star properties can be satisfactorily described 
by EOS microscopically derived from NN interactions. 
Before proceeding, it is important to also check if such 
EOS can satisfactorily describe the properties of normal 
nuclear matter. The EOS and NN interaction one uses 
for neutron star should also be applicable to normal nu- 
clear matter. As discussed in [23], many-body calcula- 
tions for normal nuclear matter using two-body NN in- 
teractions alone are generally not capable in reproducing 
empirical nuclear matter saturation properties. To rem- 
edy this shortcoming one needs to consider three-body 
forces and /or NN interactions with in-medium modifica- 
tions [23,[2l,[3|. 

A central result of the BR scaling is that the masses 
of mesons in nuclear medium arc suppressed (dropped) 
compared to those in free space [26, 27, 28, 29]. Nucleon- 
nucleon interactions are mediated by meson exchange, 
and clearly in-medium modifications of meson masses 
can significantly alter the NN interaction. These modi- 
fications could arise from the partial restoration of chi- 
ral symmetry at finite density/temperature or from tra- 
ditional many-body effects. Particularly important are 
the vector mesons, for which there is now evidence from 
both theory [H, [ij, Q and experiment [H, HI] that the 
masses may decrease by approximately 10 — 15% at nor- 
mal nuclear matter density and zero temperature. For 
densities below that of nuclear matter, a linear approx- 
imation for the in-medium mass decrease has been sug- 
gested namely 

^ = 1-C^, (11) 
my no 

where my is the vector meson mass in-medium, n is the 
local nuclear matter density and no the nuclear matter 
saturation density. C is a constant of value ^ 0.10 — 
0.15. BR scaling has been found to be very important for 
nuclear matter saturation properties in the ring -diagram 
calculation of symmetric nuclear matter of [22]. 

It is of interest that the effect of BR scaling in nu- 
clear matter can be well represented by an empirical 
Skyrme three-body force. (2^ The Skyrme force has been 
a widely used effective interaction in nuclear physics and 
it has been very successful in describing the properties 
of both finite nuclei and nuclear matter [43. It has both 



two-body and three-body terms, namely 

W™e=^n*,j)+ ^3b{l,J,k). (12) 

i<j i<j<k 

Here V{i,j) is a momentum dependent zero-range inter- 
action. Its three-body term is also a zero-range interac- 
tion 

Vsbihj, k) = h5{ri ~ rj)S{rj - nt) (13) 

which is usually expressed as a density-dependent two- 
body interaction of the form 

Vn{l,2) = 1(1 + X3Pa)t3Siri - r2)n{rav), (14) 

where Po- is the spin-exchange operator and rav = 5 (n. + 
rij)- ^3 and are parameters determined by fitting cer- 
tain experimental data. The general structure of Vskyrme 
is rather similar to the effective interactions based on ef- 
fective field theories (EFT) [s^], with V{i,j) correspond- 
ing to Viow-k and V{i,j, k) to the EFT three-body force. 
The Skyrme three-body force, however, is much simpler 
than that in EFT. 



IV. RESULTS AND DISCUSSIONS 
A. Symmetric nuclear matter and BR scaling 

When an EOS is used to calculate neutron-star prop- 
erties, it is important and perhaps necessary to first test 
if the EOS can satisfactorily describe the properties of 
symmetric nuclear matter such as its energy per particle 
Eq/A and saturation density no. In principle, only those 
EOSs which have done well in this test are suitable for 
being used in neutron star calculation. In this subsection, 
we shall calculate properties of symmetric nuclear mat- 
ter using the low-momentum ring-diagram EOS which we 
will use in our neutron star calculations, to test if it can 
meet the above test. As described in section II, we first 
calculate the Viow-k interaction for a chosen decimation 
scale A. Then we calculate the ground state energy per 
particle Eq/A using Eq.(6) {Eq = eI'''"' + AEq.) with 
the pphh ring diagrams summed to all orders. 

In the above calculation, the choice of A plays an im- 
portant role. As discussed in f22j, the tensor force is im- 
portant for nuclear saturation and therefore one should 
use a sufficiently large A so that the tensor force is not 
integrated out during the derivation of Viow~k- Since 
the main momentum components of the tensor force has 
k ^ 2fm~^, one needs to use A ~ 3/m~^ or larger. 
There is another consideration concerning the choice of 
A. The density of neutron star interior is very high, sev- 
eral times larger than no. To accommodate such high 
density, it is necessary to use sufficiently large A, sug- 
gesting a choice of A larger than ~ 3/m~^. As dis- 
cussed in [iO], a nice feature of Viow-k is its near unique- 
ness: The Viow-kS derived from various different realistic 
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NN potentials are practically identical to each other for 
A 2.1/to~^, while for larger As the resulting Viow-kS 
begin to have noticeable differences but are still similar 
to each other for A up to about 3.5/m~^. This and the 
above considerations have led us to choose A between ~ 3 
and ^ for our present study. The dependence 

of our results on the choice on A will be discussed later 
on. 




FIG. 2: Ring-diagram EOSs for symmetric nuclear matter 
with Viom-k derived from CDBonn and Bonn A potentials and 
A=3 and 3.5/m~^. 
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FIG. 3: Ring-diagram EOSs for symmetric nuclear matter 
given by Viow-k alone, Viow-k with linear (BRi) and non- 
linear {BR2) scalings, and Viow-k plus the three-body force 
(TBF) Vsb of Eq.(13). A=3.5/m"^ used for all cases. See 
text for other explanations. 

We have carried out Viow-k ring-diagram calculations 
for symmetric nuclear matter using several NN poten- 
tials (CDBonn Nijmegen [H, Argonne V18 and 
Bonn A [l^) with several values of A ranging from 3.0 
to 3.5/to~^. In Fig. 2 we present some representative 
results using CDBonn and BonnA potentials, the results 
for other potentials and As being very similar. As shown, 
the results for small densities are nearly independent of A 



within the range considered. But for larger densities, the 
results have significant variations with A and potentials, 
this being possibly a reflection of the different short-range 
repulsions contained in the potential models. A common 
feature of our results is, as displayed in the figure, that 
the calculated Eq /A and saturation Fermi momentum kp 
are all both too high compared with the empirical values 
(Eq ~-16 MeV and fc^ - 1.35/m-i or uq -0.16/m-3). 

As discussed in section III, the above situation can be 
largely improved by way of using a Vnn with BR scaling. 
In [23] authors have carried DBHF calculations for sym- 
metric matter using a BR-scaled BonnB NN potential, 
and they obtained results in good agreement with the 
empirical values, largely improved over those from the 
unsealed potential. In [22|, ring-diagram EOS calcula- 
tions for symmetric nuclear matter have been performed 
using the Nijmegen potential without and with BR scal- 
ing, the latter giving highly improved results for nuclear 
matter saturation. It should be useful if the above ef- 
fect on nuclear saturation from BR scaling also holds for 
other NN potentials. To study this, we use in the present 
work a different potential, the BonnA potential [l^, for 
investigating the effect of BR scaling on ring-diagram cal- 
culations for symmetric nuclear matter. In Fig. 3, results 
of such ring-diagram calculations for symmetric nuclear 
matter with and without BR scaling are presented. For 
the scaled calculation, the mesons (p, oj, a) of the BonnA 
potential are slightly scaled according to Eq.(ll) with the 
choice of Cp = 0.113, = 0.128 and C^r = 0.102. These 
values are chosen so that the calculated Eq /A and kp are 
in satisfactory agreement with the empirical values. The 
EOS given by the above BR-scaled potential is shown by 
the top curve of Fig. 3 (labelled as and it has 

_Eo/^=-15.3 MeV and fc^=1.33 fm~^, in good agreement 
with the empirical values. In addition, it has compression 
modulus K=225 MeV. The result using Viow~k alone is 
also shown in Fig. 3 (bottom curve). Clearly BR scaling 
is also important and helpful for the BonnA potential in 
reproducing empirical nuclear matter saturation proper- 
ties. 

We shall now discuss if the above effect of BR scaling 
can be simulated by an empirical three-body force of the 
Skyrme type. It is generally agreed that the use of two- 
body force alone can not satisfactorily describe nuclear 
saturation, certain three-body forces are needed to en- 
sure nuclear saturation [35]. There are basic similarities 
between three-body force and BR scaling. To see this, let 
us consider a meson exchanged between two interacting 
nucleons. When this meson interacts with a third spec- 
tator nucleon, this process contributes to BR scaling or 
equivalently it generates the three-body interaction. In 
[22 ] , it was already found the ring-diagram results of BR- 
scaled Viow-k derived from the Nijmegen potential can be 
well reproduced by the same calculation except the use 
of the interaction given by the sum of the Viow-k plus the 
empirical three-body force (TBF) V^^b of Eq.(13). (Note 
that Vzb is calculated using Eq.(14) with n being the local 
nuclear matter density.) Here we repeat this calculation 
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using a different potential, namely the BonnA potential. 
The strength parameter is adjusted so that the low- 
density (<'^ no) EOS given by the {Viow-k+Vsb) calcula- 
tion are in good agreement with that from the BR-scaled 
Viow-k- (We fix the parameter of Eq.(14) as zero, 
corresponding to treating the ^5*0 and '^Si channels on 
the same footing.) Results for such a calculation, with 
^3 chosen as 2000 MeV/m^ are presented as the middle 
curve of Fig. 3 (labelled as ^Viow-k plus TBF'). As shown, 
for kp lAfm~^ they agree very well with the results 



from the BR-scaled V/n 



The above (Viow-k+Vsb) cal- 



culation gives Eo/A=-lA.7 MeV and fc^ = 1.40/to"^ in 
satisfactory agreement with the BR-scaled results given 
earlier. Its compression modulus is k=140 MeV. 

It should be noticed, however, for kp >^ the 
curve for '_Bi?i scaled Viow-k rises much more rapidly 
(more repulsive) than the ^Viow-k plus TBF' one. The 
compression modulus given by them are also quite differ- 
ent, 225 versus 140 MeV. These differences may be re- 
lated to the Hnear BR-scaling adopted in Eq.(ll). This 
scaling is to be used for density less than ~ tiq. For den- 
sity significantly larger than no, such as in the interior of 
neutron star, this linear scaling is clearly not suitable. 

To our knowledge, how to scale the mesons at high den- 
sities is still an open question [2^, |23, [H, |2^, \^ . In the 
present work, we have considered two schemes for extend- 
ing the BR scaling to higher densities: One is the above 
Skyrme-type extrapolation; the other is an empirical 
modification where in the high density region a nonlin- 
ear scaling is assumed, namely m* /m = (1 — C{n/no)^) 
with B chosen empirically. The exponent i? is 1 in the 
linear BR scaling of Eq.(ll). As seen in Fig. 3, the linear 
BR-scaled EOS agrees well with the 'Viow-k plus TBF' 
EOS only in the low-density (<^ no) region, but not so 
for densities beyond. Can a different choice of B give 
better agreement for the high-density region? As seen 
in Fig. 3, to obtain such better agreements we need to 
use a scaling with weaker density dependence than BRi . 
Thus we have considered B < \, and have found that 
the EOSs with B near 1/3 have much improved agree- 
ments with the Skyrme EOS in the high density region. 
To illutrate, we have repeated the ^BRi EOS of Fig. 
3 with only one change, namely changing B from 1 to 
0.3. (The scaling parameters Cs are not changed, for 
convenience of comparison.) The new results, labelled 
as 'i3i?2', are also presented in Fig. 3. As seen, 'BR2' 
and ^Viow-k plus TBF' are nearly identical in a wide 
range of densities beyond ^ hq. This is an interesting 
result, indicating that below ^ uq the 'VJotu-fc plus TBF' 
EOS corresponds to the linear BRi scaled EOS while 
beyond ~ uq the nonlinear BR2 one. The BRi and 
BR2 EOSs have a small discontinuity (in slope) at no, 
and the above EOS with TBF is practically a continuous 
EOS with good fitting to both. As we shall discuss later, 
the above three-body force is also important and desir- 
able for neutron-star calculations involving much higher 
densities. Possible microscopic connections between the 
Skyrme three-body force and BR scalings are being fur- 



ther studied, and we hope to report our results soon in a 
separate publication. 

The ring-diagram nuclear matter EOSs using the 
^Viow-k plus TBF' interaction are in fact rather insen- 
sitive to the choice A. As discussed earlier, a suitable 
range for A is from ~ 3 to 3.5/7n. So in carrying 
out the above calculations, one first chooses a A within 
the above range. Then ^3 is determined by the require- 
ment that the low-density no) ring-diagram EOS 
given by i?i?i-scaled Viow-k is reproduced by that from 
(Viow-k+V^b)- In Fig. 4 we present some sample results 
for A= 3 and 3.5 fm~'^ with CDBonn and BonnA poten- 
tials, aU using ^3=2000 MeV/ni^. Note that this ^3 value 
is for A=3.5 frn~^ and BonnA potential; for convenience 
in comparison it is here used also for the other three cases. 
It is encouraging to see that within the above A range our 
results are remarkably stable with regard to the choice 
of both A and ^3. The four curves of Fig. 4 are nearly 
overlapping, and their {Eo/A,kp,K) values are all close 
to (-15 MeV, 1.40/m-i, 150 MeV). We have repeated the 
above calculations for the Nijmegen and Argonne V18 po- 
tentials, and have obtained highly similar results. As we 
shall discuss in the next subsection, the inclusion of V^t 
is also important in giving a satisfactory neutron-matter 
ring-diagram EOS. Calculations of neutron star proper- 
ties using the above (Viow~k + Vsb) interaction will also 
be presented there. Unless otherwise specified, we shall 
use from now on A=3.5 fm^^ for the decimation scale 
and ^3=2000 MeV/ni^ for the three-body force Vsb- 
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FIG. 4: Ring-diagram EOS for symmetric nuclear matter with 
the interaction being the sum of Viow-k and the three-body 
force of Eq.(13). Four sets of results are shown for CDBonn 
and BonnA potentials with A=3 and 3.5 fm~^. A common 
three-body force of t3=2000 MeV/m® is employed. 



B. Neutron star with neutrons only 

As a preliminary test of our ring-diagram EOS, in this 
subsection we shall consider neutron stars as composed 
of pure neutron matter only. This simplified structure is 
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convenient for us to describe our methods of calculation. 
In addition, this also enable us to check how well can 
the properties of neutron stars be described under the 
pure-neutron matter assumption. Realistic neutron stars 
have of course more complicated compositions; they have 
nuclei crust and their interior coniposed of neutrons as 
well as other elementary particles fr, ^ . We shall study 
the effects of using /3-stable and nuclei-crust EOSs in our 
neutron star calculations in the next subsection. In the 
present work we consider neutron stars at zero tempera- 
ture. 

Using the methods outlined in section II, we first cal- 
culate the ground-state energy per particle Eg /A for neu- 
tron matter. Then the energy density e, inclusive of the 
rest-mass energy, is obtained as 

e{n) ^ + m„c^) (15) 

where c is the speed of light and m„ the nucleon mass. By 
differentiating Eq/A with density, we obtain the pressure- 
density relation 

2 d{Eo/A) 

From the above two results, the EOS s{p) is obtained. It 
is the EOS e{p) which is used in the solution of the TOV 
equations. 

To accommodate the high densities in the interior 
of neutron stars, we have chosen A=3.5 fm~^ for our 
present neutron star calculation. Our ring-diagram EOS 
for neutron matter is then calculated using the in- 
teraction {Viow-k +Vzb) with the parameter 2000 
MeV/m^. Note this value was determined for symmetric 
nuclear matter, as discussed in subsection A. Is this ^3 
also appropriate for the neutron matter EOS? We shall 
address this question here. In Fig. 5 we present results 
from the above neutron matter EOS calculations for four 
interactions (CDBonn, Nijmegen, Argonne V18, BonnA). 
It is seen that the EOSs given by them are quite close to 
each other, giving a nearly unique neutron-matter EOS. 
Friedman and Pandharipande (FP) [49[ have carried out 
variational many-body calculations for neutron matter 
EOS using the two- and three-nucleon interactions V14 
and TNI respectively, their EOS results also shown in 
Fig. 5. Brown [SO] has carried out extensive studies 
of neutron matter EOS, and has found that the FP EOS 
can be reproduced by the EOS given by certain empirical 
Skyrme effective interactions (with both two- and three- 
body parts). As seen in Fig 5, our results agree with the 
FP EOS impressively well. For comparison, we present in 
Fig. 5 also the CDBonn EOS without the inclusion of V^b 
(i.e. ^3=0). It is represented by the dotted-linc, and is 
much lower than the FP EOS, particularly at high densi- 
ties. For n <~ no/2 the effect of VJjf, is rather small, and 
in this density range one may calculate the EOS using 
Viow~k alone. Clearly the inclusion of V^b with ^3=2000 
MeV/m^ is essential for attaining the above good agree- 
ment between our EOSs and the FP one. It is of interest 



that the ^3 value determined for symmetric nuclear mat- 
ter turns out to be also appropriate for neutron matter. 

In Fig. 6, our results for the EOS e{p) are presented, 
where the EOSs given by various potentials are remark- 
ably close to each other. The inclusion of Vsb is found 
to be also important here. As also shown in Fig. 6, the 
EOS given by Viow-k alone (without V^b) hes consider- 
ably higher than those with V^b- It is of interest that 
for a given pressure, the inclusion of V^b has large ef- 
fect in reducing the energy density. We have chosen to 
use A= 3.5 fm^^, and this limits the highest pressure 
PA which can be provided by our ring-diagram EOS cal- 
culation. As shown in the figure, the highest pressure 
there is about 650 MeV/Jm'^. But in neutron star cal- 
culation we need EOS at higher pressure such as 1000 
MeV//m3 (or ~ 4 X IQ-^MocVfem^). The EOS at such 
high pressure is indeed uncertain, and some model EOS 
has to be employed. In the present work we shall adopt 
a polytrope approach, namely we fit a section of the cal- 
culated EOS near the maximum-pressure end by a poly- 
trope e(p) = ap'^ and use this polytrope to determine the 
energy density for pressure beyond pA- (In our fitting the 
section is chosen as 0.8 to 1)pa-) The polytrope EOS 
has been widely and successfully used in neutron star cal- 
culations [ni, 13 ■ fs-ct we have found that our calcu- 
lated EOS, especially the section near its high-pressure 
end, can be very accurately fitted by a polytrope. In 
Table I we list the polytropes obtained from the above 
fitting for four NN interactions. It is seen that the four 
polytropes are close to each other. The exponent 7 plays 
an important role in determining the neutron-star maxi- 
mum mass. 
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FIG. 5: Ring-diagram neutron matter EOS obtained from 
four realistic NN potentials. The interaction ^Viow-k plus 
TBF' is used. The solid line with filled small circles repre- 
sents the results from the variational many-body calculation 
of Friedman-Pandharipande [i^ . The dotted line denotes the 
EOS using CDBonn- Vioui-fc only (TBF suppressed). 

In obtaining the neutron star properties, we numeri- 
cally solve the TOV equations (1) by successive integra- 
tions. In so doing, we need to have the pressure Pc at 
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FIG. 6: Neutron matter s{p) obtained from four realistic NN 
potentials. The upper-left thin line denotes the e(p) from 
CDBonn- Vjo^_fc only (TBF suppressed). 



TABLE I: Fitted polytrope ap'^ for high pressure region. See 
text for other explanations. 



Potentials 


a* 


7 


CDBonn 


69.69 ± 1.01 


0.4876 ± 0.0022 


Nijmegen 


69.99 ±1.01 


0.4885 ± 0.0021 


BonnA 


72.30 ± 1.01 


0.4779 ± 0.0021 


Argonne V18 


67.71 ± 1.01 


0.4887 ± 0.0021 



the center of the neutron star to begin the integration. 
As wc shall sec soon, different PcS will give, e.g. different 
masses for neutron stars. Wc also need the EOS e{p) for 
a wide range of pressure. As discussed above, we shall 
use the ring-diagram EOS for pressure less than pA and 
the fitted polytrope EOS for larger pressure. In Table 
II, we list some typical results for the neutron star mass 
M and its corresponding radius R and static moment of 
inertia I. (The calculation of / will be discussed later.) 
They were obtained with four different center pressures 
Pc, and as seen these properties of the neutron star vary 
significantly with Pc- 



TABLE 11: Neutron stars with different center pressures. 



PclMQ&'/km^] 


M[Mq] 


R[km] 




8.07 X 10" 


-V 


0.101 


11.58 


3.78 


7.18 X 10" 


-6 


0.347 


10.12 


13.51 


5.38 X 10" 


-5 


1.037 


10.10 


50.02 


2.33 X 10" 


-4 


1.597 


10.00 


70.69 



We present some of our calculated results for the mass- 
radius trajectories of neutron stars in Fig. 7. They were 

obtained using the CDBonn Viow-k (A=3.5 fm~^) with 
and without the three-body force Vsb (^3=2000 MeV/m'^) 
discussed earlier. As seen, the inclusion of V^b signifi- 
cantly increases both the maximum neutron star mass M 
and its corresponding radius R; the former increased from 
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FIG. 7: Mass-radius trajectories of pure neutron stars from 
ring-diagram EOSs given by the CD-Bonn Viom-k interaction 

with and without the three-body force (TBF) V'm,- Only stars 
to the right of maximum mass arc stable against gravitational 
collapse. Causality limit is indicated by the straight line in 
the upper left corner. 




FIG. 8: Density profiles of maximum- mass neutron stars of 
Fig. 7. 



~ 1.2 to ~ I.8M0 and the latter from ~ 7 to ~ 9km. The 

above results are understandable, because V36 makes the 
EOS stiffer and consequently enhances both M and R. 
Note that our results are within the causality limit. We 
have repeated the above calculations using the Nijmegen, 
Argonne and BonnA potentials, with results quite simi- 
lar to the CDBonn ones. In Fig. 8, we present the den- 
sity profiles corresponding to the maximum-mass neutron 
stars of Fig. 7. It is clearly seen that the inclusion of the 
three-body force TBF has important effect in neutron- 
star's density distribution, reducing the central density 
and enhancing the outer one. 

We have also performed calciilations using the BRi- 
scaled Viow-k interaction (BonnA and A=3.5 fm~^) 
without Vsb- The resulting maximum mass and its ra- 
dius given are respectively ^ 3.2Mq and ~ 12km, both 
being considerably larger than the values of Fig. 7. This 
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is also reasonable, because, as was shown in Fig. 3 the 
i3i?i-scaled EOS is much stiffer than the ^Viow-k plus 
TBF' one. It may be mentioned that if the neutron- 
matter EOS given by the BR-scaled interaction is plot- 
ted in Fig. 5, it would be very much higher, especially in 
the high density region, than the FP EOS shown there. 
However, the 'Viow-k plus TBF' ones are very close to 
the FP one as shown earlier. We feel that the above 
comparison is a further indication that the linear BRi 
scaling of Eq.(ll) is not suitable for high density. It is 
suitable only for density up to about ^ uq. 

Moment of inertia is an important property of neutron 
stars. [H, [111 Here we would like to calculate this quantity 
using our Viow-k ring diagram formalism. Recall that we 
have used the TOV equations (1) to calculate neutron 
star mass and radius, and in so doing we also obtain 
the density distribution inside neutron stars. From this 
distribution, the moment of inertia / of neutron stars is 
readily calculated. It may be noted that the TOV equa- 
tions are for spherical and static (non-rotating) neutron 
stars, and the I so obtained is the static one for spherical 
neutron stars. The moment of inertia for rotating stars 
are more complicated to calculate, but for low rotational 
frequencies (less than ~ 300Hz) they are rather close to 
the static ones [s^. In Fig. 8, we present our results for 
two calculations, the interactions used being the same as 
in Fig. 7. It is seen that the the inclusion of our three- 
body force V^b (TBF) largely enhances the moment of 
inertia of maximum-mass neutron star. 

The measurement of neutron-star moment of inertia is 
still rather uncertain, and the best determined value so 
far is that of the Crab pulsar (97± 38 MqUtu^) [U. For 
M > I.OMq, Lattimer and Schutz [s^j have determined 
an empirical formula relating the moment of inertia / of 
neutron stars to their mass AI and radius R, namely 

/ w (0.237 ±0.008)Afi?2 

. ^ M km M km, A. 

To check if our calculated {M, R, I) are consistent with 
this empirical relation, we have computed / using our cal- 
culated M and R values (with TBF as in the top curve 
of Fig. 9) as inputs to Eq.(17). Results of this compu- 
tation are also shown in Fig. 9. As shown, they are in 
good general agreement with the empirical formula. Es- 
pecially our moment of inertia at maximum mass agrees 
remarkably well with the corresponding empirical value. 
We have also repeated the above computation with other 
potentials and obtained similar results. 



C. Effects from /3-stable and nuclei-crust EOSs 

In the preceding subsection, we considered neutron 
star as composed of neutrons only, and we have obtained 
rather satisfactory results. Would the quality of them be 
significantly changed when we use a more realistic com- 
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FIG. 9: Pure neutron stars' moments of inertia calculated 
from CDBonn potential with and without three-body force 
(TBF). Results from the empirical formula (17) for M > 
I.OM0 with TBF are denoted by the solid line. 

position? As a small-step improvement, in this subsec- 
tion we shall first carry out calculations using the ring- 
diagram /9-stable EOS composed of neutrons, protons, 
electrons and muons only. The results of them will be 
briefiy compared with those obtained with neutrons only. 
Calculations using a combination of the BPS EOS [34] 
inside the nuclei crust and our /3-stable EOS for the in- 
terior will also be carried out. The crust of the neutron 
star is c omp osed of two parts, the outer and inner crust 
m m, [EE HE HI- The choice of the density regions 
defining these crusts and how to match the EOSs at the 
boundaries between different regions will be discussed. 

Let us first discuss our /3-stable EOS, where the com- 
position fractions of its constituents are determined by 
the chemical equilibrium equations 

^71=^? + Me (18) 

fJ-e = (19) 

together with the charge and mass conservation condi- 
tions 

rip = Tie + (20) 
n — Tin + Up + rie + n^. (21) 

In the above, fj,n, IJ^p, l^e and /x^ are the chemical po- 
tentials for neutron, proton, electron and muon respec- 
tively, and their densities are respectively n„, Up, Ue and 
n^. The total density is n. For a given n, the com- 
position fractions of these constituents are determined 
by the above equations. Note that these equations are 
solved self-consistently, since the chemical potentials and 
densities are inter-dependent. We have used iteration 
methods for this solution. Clearly the composition frac- 
tions are not uniform inside the /3-stable neutron star; 
they depend on the local density. For example, the pro- 
ton fractions x = n^/n in different density regions of /3- 
stable neutron star are generally different. In solving the 
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above equations, we have calculated the chemical poten- 
tials /i„ and /j,p using the HF approximations as indicated 
by Eq.(8). Since the interactions involving electrons and 
muons are much weaker than the strong nucleon ones, 
we have treated them as free Fermi gases and in this way 
their chemical potentials are readily obtained. 

Results for the proton fractions x calculated from the 
above equations are displayed in Fig. 10; they were ob- 
tained using four Viou)-fc(A— 3.5 fm^^) interactions, all 
with the three-body force y3f,(i3=2000 MeV/m^). We 
note that our calculated proton fractions are all quite 
small, the maximum proton fraction given by BonnA po- 
tential being ~ 7% and that by both Argonne V18 and 
Nijmegen potentials being ^ 2%. Also they exhibit a 
saturation behavior, being near maximum at density be- 
tween ~0.6 and ^0.8 fm^'^ and diminishing to near zero 
on both sides. Our results suggest that /3-stable neutron 
star has small proton admixtures only within a interme- 
diate layer; the neutron star's core and surface layer are 
both essentially pure neutron matter. Proton fractions in 
/3-stable neutron stars are an important topic and have 
been extensively studied 0, [El] . They are closely related 
to the density dependence of symmetry energy, which 
is being determined in several laboratories There 
have been a number of calculations for these fractions 
using different many-body methods and different inter- 
actions; their results are, however, widely different from 
each other (see Fig.l of ^^), some of them being close 
to ours. Further studies of the proton fractions would 
certainly be in need and of interest. 
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FIG. 10: Proton fraction of /3-stable neutron star from re- 
alistic NN potentials. Symbols are BonnA(*), CDBonn(o), 
Argonne V18 (□) and Nijemgen (x). The interaction ^Viow-k 
plus TBF' is used. 

With our calculated proton fractions, we proceed to 
calculate the properties of /3-stable neutron stars whose 

-e^, where n is the 
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energy-density EOS is e{n) — e 
total density of Eq.(21). Electrons and muons are treated 
as free Fermi gas and their energy densities are readily 
obtained. e„p is the neutron-proton energy density to 
be evaluated using the ring-diagram method described 



in section II. This energy density is in fact Snp^nnpTo) 
where n^p = Un + rip, the combined nucleon density, 
and a is the asymmetric parameter of Eq.(5). Calcu- 
lations for /9-stable neutron stars are more complicated 
than the pure neutron matter case of subsection B, for 
which a ~ 1 independent of the total density n. In con- 
trast, here we need to calculate e„p for many (n„p,Q;) 
values since they are dependent on n (see Eqs. (18-21)). 
Then the EOS e„p(p), which expresses energy density in 
terms of pressure p, is obtained by density differentia- 
tions of enpirinp, ct), similar to what we did in subsection 
A. Then by solving the TOV equations, the various prop- 
erties of /3-stable neutron stars are obtained. 
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FIG. 11: Mass-radius trajectories of neutrons stars obtained 
using only the ring-diagram /3-stable EOS (/3-alone) and a 
combination of the ring and nuclei-crust EOSs with inner- 
crust boundary nt= 0.04 (/3-crusti) and 0.05 fm~^ (/3-crust2). 
The 'CDBonn-Vlou)-fe plus TBF' potential is used for the ring 
EOS. See the caption of Fig. 7 for other explanations. 

To illustrate our results for the mass-radius trajectories 
for the /3-stable neutron stars consisted of (n,p, e,/i), we 
have performed such a calculation for the CDBonn poten- 
tial and present its results in Fig. 11 (labelled '/3-alone'). 
As seen it is quite similar to the corresponding one of Fig. 
7 for the pure neutron EOS. The trajectories using the 
same method but with other NN potentials (Nijmegen, 
Argonne and BonnA) have also been calculated, and are 
also very similar to the corresponding pure-neutron ones. 
This close similarity indicates that the effect from the 
admixture of {p, e, fj.) is not important as based on our 
present calculation, this being largely due to the small- 
ness of the proton fractions discussed earlier. 

So far we have carried out microscopic calculations of 
neutron stars with the assumption that they are made 
of a homogeneous medium composed of neutrons, pro- 
tons, electrons and muons. It is believed that the crust 
of neutron stars is not of such homogeneous medium; it 
is " nuclei crust" where nucleons are clustering into nuclei 
[3^ . [s^ l . Here we would like to make some estimates on 
the nuclei-crust corrections to our present calculations. 
In our estimates, we employ three different EOSs for 
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TABLE III: Maximum mass and the corresponding radius 
and moment of inertia of /3-stable neutron stars with nuclei- 
crust boundary nt=0.04 fm~^. The three-body force V^b is 
included for the results in the first four rows, but is not in the 
last. 



Potentials 


M[Mq] 


R[km] 


I[Mcs,km^] 


CDBonn 


1.80 


8.94 


60.51 


Nijmegen 


1.76 


8.92 


57.84 


BonnA 


1.81 


8.86 


61.09 


Argonne V18 


1.82 


9.10 


62.10 


CDBonn(T/36=0) 


1.24 


7.26 


24.30 
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FIG. 12: Mass-radius trajectories of neutron stars calculated 
with a combination of /3-stable ring-diagram EOS for the core 
and the nuclei-crust EOSs for the crusts. Ring-diagram EOSs 
given by four NN potentials, all with the Vsb three-body force, 
A=3.5/m"^ and t3=2000 MeVfm^ are used. nt=0.04 /m"^ 
is used for the inner crust boundary. See the caption of Fig. 
7 for other explanations. 
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FIG. 13: Moments of inertia of neutron stars of Fig. 12. 



the outer-, inner-crust and core regions. These regions 
refer to the density regions n < Uout, "-out < n < nt 
and n > nt respectively, with Uout = 2.57 x lQ~^fm~^ 
[sl . [56j . nt is the transition density separating the in- 
ner crust and the homogeneous core, and several models 
have been employed to determine its value [13, HI] . For 
the outer-crust region we use the well-known BPS nuclei- 
EOS (3^. For the core region, our /3-stable ring-diagram 
EOS will be employed. The EOS in the inner-crust region 
is somewhat uncertain, and so is the transition density 
separating the inner crust and core. We shall use in our 
calculations nt — 0.05 and 0.04/m"'^ [54, 55J to illus- 
trate the effect of the nuclei crust. Following 0,(53], we 
use in the inner-crust region a polytropic EOS, namely 
p ~ a + be^^^ with the constants a and b determined by 
requiring a continuous matching of the three EOSs at 
nout and nt. 



In Fig. 11, our results for the mass-radius trajecto- 
ries using the above three EOSs with rit=0.04 and 0.05 
/to~^, labelled /3-crusti and /3-crust2 respectfully, are 
compared with the trajectory given by the /3-stable alone 
(namely nt—0). As seen, the effect from the nuclei-crust 
EOS on the maximum neutron-star mass and its radius 
is rather small, merely increasing the maximum mass by 
^ 0.02Mq and its radius ~ 0.1km as compared with 
the P-alone ones. However, its effect is important in the 
low-mass large-radius region, significantly enhancing the 
neutron-star mass there. That the maximum mass is not 
significantly changed by the inclusion nuclei-crust EOS 
is consistent with Fig. 8 which indicates the mass of 
maximum-mass neutron stars being confined predomi- 
nantly in the core region. It may be mentioned that our 
ring-diagram EOS is microscopically calculated from re- 
alistic NN interactions, while the crust EOSs are not. 
So there are disparities between them. It would be use- 
ful and of much interest if the crust EOSs can also be 
derived from realistic NN interactions using similar mi- 
croscopic methods. Further studies in this direction are 
needed. 

In. Fig. 12 we present our mass-radius results using the 
above three EOSs with nt=0.04 fm~^. Four NN poten- 
tials are employed, and they give similar trajectories, es- 
pecially in the high- and low-mass regions. A correspond- 
ing comparison for the moment of inertia is presented in 
Fig. 13; again the results from the four potentials are 
similar. In Table III, our results for the maximum neu- 
tron mass and its radius and moment of inertia using the 
above combined EOSs arc presented, and as seen the re- 
sults for the maximum-mass neutron star given by the 
four potentials are indeed close to each other. It is also 
seen that the effect of the three-body force is quite im- 
portant for M, R and /, as illustrated by the CDBonn 
case. 



V. SUMMARY AND CONCLUSION 

We have performed neutron-star calculations based on 
three types of EOSs: the pure-neutron ring-diagram, 
the /3-stable (n,p, e, /i) ring-diagram and the BPS nuclei- 
crust EOS. The ring-diagram EOSs, where the pphh ring 
diagrams are summed to all orders, are microscopically 
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derived using the low-momentum interaction Viow-k ob- 
tained from four realistic NN potentials (CDbonn, Ni- 
jmegen, Argonne V18, BonnA). We require that the EOS 
used for neutron stars should give satisfactory saturation 
properties for symmetric nuclear matter, but this require- 
ment is not met by our calculations using the above po- 
tentials as they are. Satisfactory nuclear matter satura- 
tion properties can be attained by using the above po- 
tentials with the commonly used linear BR scaling (BRi) 
where the the masses of in-medium mesons are slightly 
suppressed compared with their masses in vacuum. But 
this linear scaling is not suitable for neutron stars; the 
maximum mass of neutron star given by our BRi ring- 
diagram calculation is ^ 3.2Mq which is not satisfactory. 
BRi is suitable only for low densities; it needs some ex- 
tension so that it can be applied to the high densities 
inside the neutron star. We have used an extrapolation 
method for this extension, namely we add an empirical 
Skyrme-type three-body force V^b to Viow-k- We have 
found that the EOS given by this extrapolation agrees 
well with the EOS obtained from linear BRi scaling for 
low densities, but for high densities it agrees well with 
that from a nonlinear BR2 scaling. The EOS using the 
above extrapolation gives satisfactory saturation proper- 
ties for symmetric nuclear matter, and for neutron matter 
it agrees well with the FP EOS for neutron matter. 

The effects from V^b have been found to be both im- 
portant and desirable. Compared with the results given 
by the unsealed Viow-k, it increases the maximum mass 
of the neutron star and its radius and moment of iner- 
tia by - 40% , ~ 20% and ~ 150% respectively. The 
proton fractions are found to be generally small (< 7%), 



making our neutron-star results using the pure-neutron 
EOS and those using the /3-stable EOS being nearly the 
same. We have estimated the effect from the nuclei-crust 
EOSs by using a combination of three EOSs: the BPS 
EOS for the outer crust, a fitted polytropic EOS for the 
inner crust and our /3-stable ring-diagram EOS for the 
core region. The effect from the nuclei-crust EOSs on 
the maximum neutron-star mass and its radius is found 
to be rather small, as compared with those given by the 
calculation where the /3-stable EOS is used throughout. 
However, its effect is important in the low-mass large- 
radius region, significantly enhancing the neutron-star 
mass there. Using the above combined three EOSs, our 
results for neutron star's maximum mass and its radius 
and moment of inertia are, respectively, ~ I.8M0, ~ 9km 
and ^ GOMQkrn^ , all in good agreement with accepted 
values. 

How to extend the BR scaling to high densities is still 
an open question. Although we have obtained satisfac- 
tory results by using a nonlinear scaling for the high- 
density region, or equivalently a Skyrme-type three-body 
force, for the extension, it would still be certainly useful 
and interesting to explore other ways for doing so. Fur- 
ther studies in this direction would be very helpful in 
determining the medium dependent modifications to the 
NN potentials in the high-density region. 
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